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It is shown that optical pulses with an average position accuracy beyond the standard quantum limit can be 
produced by adiabatically expanding an optical vector soliton followed by classical dispersion management. 
The proposed scheme is also capable of entangling positions of optical pulses and can potentially be used for 
general continuous-variable quantum-information processing. 
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If an optical pulse consists of N independent photons, then 
the uncertainty in the pulse-center position is the pulse width 
divided by y/N, the so-called standard quantum limit HI] . The 
ultimate limit permissible by quantum mechanics, however, 
is determined by the Heisenberg uncertainty principle and is 
smaller than the standard quantum limit by another factor of 
y/N, resulting in a quantum-enhanced accuracy useful for po- 
sitioning and clock synchronization applications [2]. To do 
better than the standard quantum limit, a multiphoton state 
with positive frequency correlations and, equivalently, nega- 
tive time correlations is needed J2]. Consequently, significant 
theoretical and experimental |5|] efforts have been made 
to create such a nonclassical multiphoton state. All previous 
efforts were based on the phenomenon of spontaneous photon 
pair generation in parametric processes, limiting N to 2 only. 
The resultant enhancement can only be regarded as a proof of 
concept and is too small to be useful, considering that a large 
number of uncorrected photons can easily be obtained, with a 
standard quantum limit orders of magnitude lower than the ul- 
timate limit achievable by two photons. It is hence much more 
desirable in practice to be able to enhance the position accu- 
racy of a large number of photons. In this Letter, for the first 
time to the author's knowledge, a scheme that produces a mul- 
tiphoton state with positive frequency correlations among an 
arbitrary number of photons is proposed, thus enabling quan- 
tum position accuracy enhancement for macroscopic pulses as 
well. The scheme set forth therefore represents a major step 
forward towards the use of quantum enhancement in future 
positioning and clock synchronization applications. 

The proposed scheme exploits the quantum properties of 
a vector soliton, in which photons in different optical modes 
are bound together by the combined effects of group-velocity 
dispersion, self-phase modulation, and cross-phase modula- 
tion [g]. A quantum analysis shows that the average position 
of the photons in a vector soliton is insensitive to the optical 
nonlinearities and only subject to quantum dispersive spread- 
ing, while the separations among the photons is controlled 
by the balance between dispersion and nonlinearities. These 
properties are in fact very similar to those of scalar solitons 
O, [8|], so the idea of adiabatically compressing scalar soli- 
tons for momentum squeezing [9] can be similarly applied to 
vector solitons. To produce negative time correlations, how- 
ever, adiabatic soliton expansion should be performed instead. 
Given the past success of experiments on scalar quantum soli- 



tons iflOri and vector solitons 1 11], the scheme set forth should 
be realizable with current technology. The formalism should 
apply to spatial vector solitons as well, so that the position ac- 
curacy of an optical beam can be enhanced [12]. Moreover, 
the proposed scheme is capable of creating temporal Einstein- 
Podolsky-Rosen (EPR) entanglement II 1311 among the pulses 
in a vector soliton, suggesting that the vector soliton effect, 
together with quantum temporal imaging techniques J3], may 
be used for general continuous-variable quantum-information 
processing 11411 . 

For simplicity, only vector solitons with two optical modes, 
such as optical fiber solitons with two polarizations, are an- 
alyzed in this Letter, although the results can be naturally 
extended to multimode vector solitons, such as those stud- 
ied in Refs. B15I1 . Two-mode vector solitons are classically 
described by the coupled nonlinear Schrodinger equations 
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2c{\V\ 2 +fi|t/| 2 )V, where U and V are complex envelopes of 
the two polarizations, assumed to have identical group veloc- 
ities and group-velocity dispersion, t is the propagation time, 
Z is the longitudinal position coordinate in the moving frame 
of the pulses, b is the group-velocity dispersion coefficient, c 
is the self-phase modulation coefficient, and Be is the cross- 
phase modulation coefficient. For example, B = 2/3 for lin- 
ear polarizations in a linearly birefringent fiber 11611 . B = 2 
for circular polarizations in an isotropic fiber [17], and B = 1 
describes Manakov solitons lfl8ll . realizable in an elliptically 
birefringent fiber II 1611 . be < is required for solitons to exist. 
The coupled nonlinear Schrodinger equations can be quan- 
tized using the Hamiltonian H = h f dz [b( ^ f| + ^ ^ ) + 
c(U f U^UU + y + y + y1> + 2BU Jf VWV)), where U and V are 
photon annihilation operators of the two polarizations and the 
daggers denote the corresponding creation operators. The 
Heisenberg equations of motion derived from this Hamilto- 
nian are analyzed using perturbative techniques by Rand et al. 
lfl9ll . who study the specific case of Manakov solitons, and by 
Lantz et al. H20fl and Lee et al. 112111 . who numerically inves- 
tigate the photon number entanglement in higher-order vector 
solitons. As opposed to these previous studies, in this Let- 
ter the exact quantum vector soliton solution is derived in the 
Schrodinger picture, in the spirit of the scalar soliton analyses 
in Refs. Ella. 

Since the Hamiltonian conserves photon number in each 
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mode and the average momentum, one can construct si- 
multaneous Fock and momentum eigenstates with the Bethe 
ansatz \n,m,p) = -j=-fd n xd m yf nmp (xi,...,x n ,yi,...,y m ) x 

J/ 1 " (xi ) . . .J7 f (x„ )V f (y i ) . . . V f (y m ) |0) QUI, where n and m are 



the photon numbers in the two polarizations and p is the av- 
erage momentum. Using the Schrodinger equation E\*¥) = 
H |\P), one obtains 
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The soliton solution of Eq. (Q]) is 
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where C nm is a normalization constant. The energy can be calculated by substituting Eq. (01 into Eq. (HJ and is given by 
Enmp = ftbNp 2 — [n(n 2 — 1) + m(m 2 — 1 ) + 3B 2 nm(n + m)], where N = n + m. A physical state should contain a distribution 
of momentum states, say, a Gaussian, such that the time-dependent multiphoton probability amplitude is now given by 
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where Ap is determined by initial conditions and a constant 
energy term that does not affect the position and momentum 
properties of a Fock state is omitted. Although a more real- 
istic soliton state should have a superposition of Fock states 
resembling a coherent state [7], the Fock components of a co- 
herent state for A >> 1 have photon numbers very close to the 
mean value, so a Fock state should be able to adequately repre- 
sent the position and momentum properties of a coherent-state 
soliton. 

The multiphoton amplitude f nm consists of two compo- 
nents; a dispersive pulse-center component given by Eq. (0]i 
that governs the quantum dispersion of the average photon po- 
sition jji'E.jXj+llkyk)' and a bound-state component given by 
Eq. (0 that fixes the distances among the photons via the at- 
tractive Kerr potentials. It follows that the momentum-space 
probability amplitude, defined as the A-dimensional Fourier 
transform of /„„,, also consists of an average momentum com- 
ponent and a bound-state component that governs the relative 
momenta among the photons. 

If one increases b or reduces c adiabatically, the multipho- 
ton amplitude would remain in the same form, but with in- 
creased uncertainties in the relative distances as well as re- 
duced uncertainties in the relative momenta. More crucially, 
the average momentum uncertainty remains unaffected, lead- 



ing to a multiphoton state with positive momentum corre- 
lations. The adiabatic approximation remains valid if the 
change happens over a propagation time scale T » h/\E(t — 
T) —E(t = 0)|, which is on the order of the initial soliton 
period divided by N. As optical fiber solitons can typically 
propagate for a few soliton periods before loss becomes a fac- 
tor, the desired adiabatic expansion should be realizable with 
current technology. In the following it is assumed for sim- 
plicity that only c is adiabatically varied. Mathematically, in 
the limit of vanishing c, the bound-state component becomes 
relatively flat, and f nm becomes solely governed by the pulse- 
center component, 
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In the momentum space, as the bandwidth of the relative mo- 
menta is reduced and becomes much smaller than the band- 
width of the average momentum, the wavefunction in terms 
of momentum eigenstates becomes 
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where \n p ,m p ) denotes a momentum eigenstate with momen- 
tum p and n and m photons in the respective polarizations. 
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Except for the dispersive phase term, Eq. dT) is precisely the 
desired coincident-frequency state that can achieve the ulti- 
mate limit of average position accuracy ||2|], as frequency is 
trivially related to momentum via the dispersion relation. If 
the pulse is sent across one channel only, adiabatic control of a 
scalar soliton would already suffice for the purpose of tempo- 
ral uncertainty reduction, but the use of a vector soliton allows 
quantum-enhanced pulses to be sent across different channels, 
as originally envisioned by Giovannetti et al. 0], for addi- 
tional security. The same operation of position squeezing on a 
scalar soliton is previously considered by Fini and Hagelstein, 
who nonetheless dismiss this possibility due to the detrimental 
effect of quantum dispersion [9]. 

Fortunately, quantum dispersion, like classical dispersion, 
can be compensated with classical dispersion management. If 
the vector soliton propagates in another linear waveguide with 
an opposite group-velocity dispersion b', such that bt = —b't', 
where t is the propagation time in the first waveguide and t' 
is the propagation time in the second waveguide, then the dis- 
persive phase term in Eq. (0), —ibNp 2 t, can be cancelled, thus 
restoring the minimum uncertainty in the average photon po- 
sition, while the pulse bandwidth remains constant because 
the second waveguide is linear. The complete proposed setup 
is sketched in Fig.Q] To apply the scheme set forth to a spatial 
vector soliton, negative refraction 1E3I1 is required to compen- 
sate for the quantum diffraction instead. 

Pulses with 
Coincident Frequencies 



Adiabatically Increasing Dispersion or 
Decreasing Nonlinearity 




-» | Nonlinear Waveguide Dispersion Compensatt 



FIG. 1: Proposed setup of generating multiphoton states with 
quantum-enhanced average position accuracy. 

In order to understand how the quantum vector soliton so- 
lution corresponds to a classical soliton in typical experi- 
ments, consider the specific case of a Manakov soliton, where 
B = 1 . Other vector solitons should have very similar prop- 
erties given the similarity of the solutions. If the photon 
position variables are re-indexed in the following new nota- 
tions {zi,—,Zn} — {jci, ... ,x„,yi,... ,y m }, the multiphoton am- 
plitude in Eqs. and (0 becomes 
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Intriguingly, this solution is exactly the same as the scalar 
soliton solution 0], or in other words, a Manakov soliton is 
quantum-mechanically equivalent to a scalar soliton. This 



equivalence explains the discovery by Rand et al. that the 
squeezing effect of a Manakov soliton has the same opti- 
mum as a scalar soliton lfl9ll . Moreover, C nm can now be 
borrowed from the scalar soliton analysis and is given by 
C, m = \{N- \)\\c/b\ N - l /{2n)} l l 2 0]. The knowledge of C nm 
allows one to calculate the correlations among the photon po- 
sitions using standard statistical mechanics techniques. An 
expression for (Ej<; \zj — Zi\) can be derived, and, by sym- 
metry, 



(\zj-Zi\) = 
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As expected, the mean absolute distance between any two 
photons is on the order of the classical soliton pulse width, 
Wo ~ \2b/(Nc)\ 0]. Next, assume that the variance of the 
relative distance is related to the square of the mean absolute 
distance by a parameter q, 
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While an explicit expression for q is hard to derive, q must 
depend only on N by dimensional analysis, must be larger 

than 1 because (\zj — Zi\ 2 ) > (\zj — Zi\) , and is likely to be 
on the order of unity, as will be shown later. By symme- 

try, (\ Zj -z,\ 2 ) = (z 2 )-2(z l z J ) + (z 2 )=2(z 2 )-2( Zi z J ). 
Eq. ([Tol l then gives 
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Furthermore, the variance of £ ; Zj is simply given by 
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where Az 2 = l/(4Ap 2 ) +4(bNApt ) 2 . From Eqs. (TTTb and (fT2l 
the covariances can be obtained explicitly, 
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A quantum soliton solution best resembles a classical initial 
condition with independent photons when the initial covari- 
ance is zero, <Z/Z./) ;=0 = l/{4N 2 Ap 2 ) -2qb 2 /(N 3 c 2 )\ t=Q = 0, 
and 
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Incidentally, the average momentum uncertainty Ap is at the 
shot-noise level when the photons are initially uncorrelated. 
This justifies the assumption that q is on the order of unity. An 
initial condition with independent photons would then mostly 
couple to a soliton state with Ap given by Eq. ( fT4b . while cou- 
pling to continuum states should be negligible. Adiabatically 



4 



increasing \b/c\ thus makes (ziZj) negative and therefore in- 
troduces the necessary negative correlations among the pho- 
ton positions. To investigate the magnitude of the quantum 
enhancement in practice, it is useful to compare the quan- 
tum theory to the classical soliton theory, as the two regimes 
should converge when N » I and y « \/~N. According to 
the classical theory, if the ratio between the final and initial 
values of \b/c\ is y, the pulse bandwidth is also reduced by a 
factor of 7, from ~ 1 /Wo to ~ 1/ (yWb). Since the final av- 
erage position uncertainty is the same as the input value, the 
accuracy enhancement over the standard quantum limit, for 
the same reduced bandwidth ~ 1/ (yWb), is hence also given 
by y, in the regime of moderate pulse expansion y < < \/N. 
Because the ultimate soliton state, given by Eqs. d6) and (O, 
has a bandwidth given by Eq. (TBI . ~ 1 /(VNWo), the ultimate 
limit is reached only when y >> V^V- 

As the photons across different optical modes become cor- 
related via the cross-phase modulation effect, entanglement is 
expected among the pulse positions in a vector soliton. To es- 
timate the magnitude of the entanglement in terms of macro- 
scopic position variables, consider again the case of Man- 
akov solitons. Let the pulse-center coordinates of the respec- 
tive polarizations be X and Y, defined as X = ^L" =1 Z; and 

^ = Y*k=m+i z k- If n = m = N/2 is assumed for simplic- 
ity, the following statistics for X and Y can be derived using 
Eqs. O, 

Similar to a two-photon vector soliton fl], the average po- 
sition of the two pulses is affected by quantum dispersion, 
while the relative distance is bounded by the Kerr effect. 
For two initially uncorrelated pulses, the two expressions in 
Eq. $15[ have the same value. If, however, b and c are adia- 
batically manipulated, then the nonlocal uncertainty product 
((X - Y) 2 ) ((P x - Py) 2 ), where P x and P Y are the conjugate 
momenta, can remain constant under the adiabatic approxi- 
mation, while ((X-Y) 2 ) and ((P x ~Py) 2 ) can be arbitrar- 
ily varied. Since ((Px+Py) 2 ) always remains constant and 
((X + F) 2 ) can also remain the same as the input value if 
quantum dispersion is compensated, ((X — F) 2 ) ( (Px + Py) 2 ) 
or ( (X + F) 2 ) / (Px — Py) 2 ) can be arbitrarily reduced, thus re- 
sulting in EPR entanglement. Combined with quantum tem- 
poral imaging techniques, which are able to temporally re- 
verse, compress, and expand photons in each mode 0], adi- 
abatic vector soliton control potentially provides a powerful 
way of fiber-based continuous-variable quantum-information 
processing 11411 . 
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ERRATUM 

On closer inspection, it is found that the soliton solution 
given by Eq. (2) and the expression for the energy E mnp below 
Eq. (2) in the published Letter (H are valid only for the Man- 
akov soliton (B = 1) or uncoupled scalar solitons (B = 0). This 
can be seen by observing that the energies calculated using 
Eqs. (1) and (2) in Ref. [ 1] in different regions of the configu- 
ration space are different unless B = 1 or B = 0. For example, 
for n = 2 and m=l, the energy in the regions of x\ <yi <X2 
and X2 < yi < x\ calculated using Eq. (2) is different from the 
energy in the other regions, unless B — 1 or B = 0. 

This error does not affect the rest of Ref. that focuses 
on the Manakov soliton. Furthermore, it can be argued on 
physical grounds that the proposed scheme should still work 
for any positive B, despite the lack of an explicit soliton so- 
lution in general. First, it can be shown using the "center- 
of-mass" coordinate system [2] that all stable solutions of the 
Schrodginer equation given by Eq. ( 1 ) in Ref. [ 1 ] must consist 
of a dispersive center-of-mass component and a bound-state 
component. For any positive B, the bound-state component 
must become flat in the ultimate limit of adiabatic pulse ex- 
pansion, so any stable solution still approaches the form of 
Eqs. (6) and (7) in Ref. yj]. Second, the average timing jit- 
ter of the two pulses is affected only by dispersion, regardless 
of B. If the net dispersion is zero and the fiber is assumed 
to be lossless, the output jitter is the same as the input jit- 



ter. As long as the solitary wave propagation is stable against 
the adiabatic change in fiber parameters, the bandwidth can 
be adiabatically reduced, leading to a raised timing standard 
quantum limit. Hence the final jitter must be lower than the 
raised standard quantum limit regardless of B. Finally, Eq. (1) 
in Ref. [1] shows that the cross-phase modulation effect man- 
ifests itself as an attractive potential among photons across 
the two polarizations. In a stable solution of Eq. (1), the rel- 
ative position and relative momentum of the two pulses can 
therefore be controlled by adjusting the balance between dis- 
persion and cross-phase modulation. Since temporal entan- 
glement has been proved to occur for the Manakov soliton in 
Ref. this physical picture should remain qualitatively the 
same for a different cross-phase modulation coefficient, and 
temporal entanglement should also occur for any positive B. 

There are also two typographical errors in [1]. The fac- 
tor N ^_^ in Eq. (9) should read N ^_^ , and in the second- 
to-last paragraph, Y — (1/m) Ylk=m+\ z k should read Y — 
(1/w) Y!k=„+\ Zk- These typographical errors do not affect the 
rest of the Letter in any way. 
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